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Abstract 

We derive a semiclassical formula for the tunneling current of electrons 
trapped in a potential well which can tunnel into and across a wide quan- 
tum well. The calculations idealize an experimental situation where a strong 
magnetic field tilted with respect to an electric field is used. The result- 
ing semiclassical expression is written as the sum over special periodic orbits 
which hit both walls of the quantum well and are perpendicular to the first 
wall. 

PACS numbers: 03.65.Sq, 05.45.+b 
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The connection between quantum properties and the underlying classical mechanics has 
attracted wide attention in the last years (see e.g. flj and references therein). The central 
result of this investigation, the Gutzwiller trace formula, relates long-range fluctuations in 
the density of quantized levels with classical periodic orbits. Each periodic orbit of the 
classical dynamical system with period T p generates regular maxima in the level density 
separated by an energy AE P = h/T p . 

The influence of classical motion on quantum mechanical systems has also been inves- 
tigated experimentally [0- [[|. In all cases the observed oscillations were attributed to 
fluctuations due to special periodic orbits. 

In this paper we discuss a particular system to which much experimental work has been 
devoted in recent years [gj - [§]. The double well potential consisting of GaAs/(AlGa)As 
resonant tunneling diode containing a wide quantum well (QW) has been used to explore a 
relationship between classical and quantum electron dynamics. In the presence of a strong 
uniform magnetic field (B) tilted with respect to the z-axis by an angle 9 and a uniform 
electric field (e) parallel to the z-axis, the system exhibits chaotic motion for certain initial 
conditions in phase space []5| - and the experiments and numerical simulations reveal 
oscillations in the tunneling current associated with certain classical periodic orbits [^| - [12|. 
Even the scarring of wave functions for this problem has been found ||. 

The main purpose of this paper is to obtain a quantitative description of the link between 
the tunneling current and periodic orbits. We shall show that in the strict semiclassical 
approximation the main contribution will be from special periodic orbits which (i) hit the 
two walls and (ii) are perpendicular to the first barrier. Unlike the Gutzwiller trace formula 
the result is not a canonical invariant. 

The Hamiltonian (in atomic units) for motion of an electron inside the system, is effec- 
tively 2-dimensional and can be written as: 

H(y,z) = ^ + ^(y cosd ~ zsind f ~ ez + U ( z )> C 1 ) 

where U(z) is a step-wise potential representing allowed - forbidden layers along z and 
m = 0.067 is the band edge mass of an electron in GaAs. 

In an over-simplified approximation (cf. |T^|) the wave function just after the barrier 



can be written as a factorized expression corresponding to an approximate separation of 
variables: 

Mv,*) = -r= ex P(r f k(z')dz')My)- (2) 
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Here k(z) is ^-momentum and ipoiu) is the lowest Landau level eigenfunction for motion in 
the magnetic field 

Mv) = (^V ex P (-^os%-y) 2 ) , (3) 

where y =< z > ta.n6 is a small shift due to the diagmagnetic term in the perturbation 
expansion. \C\ 2 = mTi and Ti is the imaginary part of the energy of the quasi-bound state 
in the first well computed without the second barrier. 
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Using the Stokes theorem one can show that the wave function inside the QW is expressed 
through its boundary values as follows 

^(y, z) = -L / (G(y, z- y', 0)dMv', 0) 



2m 

-My'^)d z G(y^;y'^))dy\ (4) 

where G(y, z\ y', z') is the Green function for the QW. One can also prove that this result 
coincides with the first order correction for wave functions in the perturbation theory on the 
tunneling amplitude similar to one used by Bardeen in [T3I . 



Assuming that the QW is sufficiently wide we shall use the standard semiclassical approx- 
imation according to which the Green function is represented by a sum over contributions 
from classical trajectories (J) connecting final x = (y, z) and initial x' = (y', z') points 
G(y, z; y', z') = J2j Gj(x, x'), and 



r.(- 2\ m 

\JT j ^ X ■ X J 



y/2n ik 2 ki 



d*S 3 



x l 2 a- 

eH'^i. (5) 



Here Sj = Sj(x,x') is the classical action calculated along the classical trajectory (j) con- 
necting initial and final points, fij is the Maslov index for this trajectory, k\ and k 2 are 
the modulus of momentum in the initial and final points respectively, and t\ and t 2 are 
coordinates perpendicular to the trajectory in these points. 

In our case, two modifications have to be made in this standard semiclassical formula, 
(i) In order to to describe the motion inside a well with finite width walls it is necessary to 
multiply the expression above by the reflection coefficient, (ii) To take into account various 
inelastic processes we shall multiply Eq. (|5]) by a damping factor exp(— ToTj) where Tj is 
the period of the trajectory and 1/Tq plays the role of a mean free time for motion in the 
QW. In the following we shall ignore the difference between different kinds of mean free 
paths and in particular the difference between the mean value of the Green function and the 
mean value of the product of the two Green functions which is important in describing real 
experiments (see e.g. [|J). 

The calculation of the Green function when one of its arguments is very close to the 
boundary of the QW requires special attention. In this case there are 2 different contributions 
corresponding to 2 trajectories with opposite value of z-momentum. The one which hits 
the boundary has to be multiplied by the reflection coefficient. Therefore G(x,x')\ z >=q = 
(1 + ri)Gj(x, xq ), and d z >G(x, x' )\ z '=o = —ip z (^ — r i)Gj(x, xq ), where n is the reflection 
coefficient from the left wall and GA Xj Xq ) is the contribution of the trajectory starting from 
the point x = (y', 0) with p z > and ending at the point x. 

The explicit expression for the Green function and the value of the wave function in 
the vicinity of the LH barrier permits one to make a semiclassical computation of the wave 
function at any point inside the QW. We get 

1>(V, z) = £ / WjWGjiy, z- y', 0)dy'. (6) 



with Wj(y') = iCv j ^p z (j)ip (y')/2m, and Vj = (1 - rx)yjpjj)jk z + (1 + n)^k z /p z (j). Here 
p z (j) is the initial ^-momentum of trajectory j and k z = k(0) with k(z) from (Q). Under a 
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natural condition \p(J) — k z \ <C \p(j)\ this expression equals 2 and we think that this value 
is a good zeroth order approximation for this quantity. 
Without the simplifying assumption (0) 



CMy')/\/k z = My', o) and k z = d z logMv', °) 

are just the boundary values of the initial wave function on the LH barrier to be determined 
either numerically or by more refined methods of multidimensional tunneling. 

Knowledge of the wave function (H) permits us to compute the current ji(x*f) = ji(yf, d) 
at the second interface of the QW (i.e. at the RH wall): 

?*(*/) = ^-(r(x f )V^(x f ) - ^(f,)V^(2,)). (7) 

When the electron hits the second barrier it has a probability of tunneling through the 
wall. In the same approximation as above the current after the barrier, (jf), differs from 
the current before it, (jj), by a transmission coefficient through this barrier (i 2 ) jf = |^2| 2 jV 
The total imaginary part, T, of the energy of a quasi-bound state in the emitter well equals 
the total current after the second barrier 

da 2 ji \t 2 \ 2 . (8) 

Substituting here the expression (|6]) for the wave function one gets an expression for T as a 
triple integral over two initial positions and one final y coordinate: 

r = H J d y d y' w jk(y,y') J G j (y f ,y)G k (y f ,y')dy f \t 2 \ 2 , (9) 

where Wj k (y,y') = (Pz(j) + P^P (k))vjj(y)wk(y') /2m and Gj(yf,y) is the contribution (fj[) of 
a trajectory (J) which starts at point (y, 0) and ends at point (yf, d). 

In the semiclassical approximation Gj is proportional to exp(iSj/h) and in the formal 
limit H — > it is natural to perform the integration over all variables by the saddle point 
method. Assuming that the boundary function ipo(y) and other quantities are smooth in 
the scale of noticeable changes of the Green function, one concludes that in such an approx- 
imation the dominant contribution to the above integrals will be given by trajectories in the 
vicinity of saddle-points trajectories for which the following three conditions are fulfilled: 

dS 3 {y fl y)/dy f - dS k (y f ,y')/dy f = 0, (10) 



dS j (y f ,y)/dy = 0, dS k (y f ,y')/dy' = 0. (11) 

The first equation means that saddle point trajectories (labeled here by j and k) should have 
the same y component of momentum: p\j\j) = p^'(k). The equality of the p z momenta 
for these trajectories then follows from energy conservation. But two classical trajectories 
passing through the same point (yf) and having the same momenta at this point can be either 
(i) exactly the same trajectory or (ii) two different paths on the same classical trajectory. 

The second pair of saddle point equations (|TT|) signifies that the saddle point trajectories 
should have zero y component of the momentum in both points at the LH wall (i.e. they 
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have to be perpendicular to the plane z = 0). The combination of these conditions leads 
to the important conclusion that in the strict semiclassical limit the tunneling current or 
probability of decay is divided into two distinct contributions: T = T( Wey i) + r( osc ). The first 
term corresponds to the interference of an arbitrary trajectory perpendicular to the plane 
z = with itself. It has no quick dependence on external fields and in real experiments it is 
effectively washed out by taking the second derivative of the current. 

The second term, which we denote by rV osc ), includes the contributions from different 
paths on the same self-retracing periodic trajectory which connect the LH and RH walls 
and is perpendicular to the plane z — 0. The classical trajectories which give the dominant 
contribution to the resonant tunneling are built from a segment of a periodic orbit which 
hits both walls plus an arbitrary number of loops (n) around this orbit. 

Expanding the actions up to the quadratic terms, expressing the derivatives of actions 
through the monodromy matrix, performing the sum over all possible repetition numbers 
one gets 



t { r\v)\ 



2 



r(osc) — ? r„- , „ 77. x 
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fl( r )i/2 e MiS( r ) - i^p{ r ) + ij) + cc, (12) 

where D(r) = m 2 i(r) and S(r) = rS p . The sum here is taken over all primitive periodic 
orbits with p y (i) = and over up to 2 points of reflection of this trajectory with the first 
wall: yo and y' Q . S p is the classical action around the periodic orbit labeled by p [jTa] , rriij (r) is 
the monodromy matrix element computed along the first barrier for r repetitions of the 
orbit considered, R p = ri to t){p) exp(— T T P ) where r^ tot ^(p) is the total reflection coefficient 

from the walls for one loop around the given orbit, |4 <0 ^(p)| 2 = \h{p, 0l 2ex P( — 2r T (/)) 
is the total tunneling probability through the second barrier weighted by damping factors, 

and T p = (7r/8) 1/2 f p ^ri^o (2/0)^0(2/0) • 

The formula, Eq. (|T2"D, is the main result of the paper. It expresses the semiclassical 
limit of the tunneling current (which is proportional to T) as a sum over special periodic 
trajectories which (i) connect the two walls of the QW and (ii) are perpendicular to the first 
barrier. 

The important ingredient in deriving the above formula was the assumption related to the 
applicability of the semiclassical approximation. Whereas for the integration over the final 
coordinate it can easily be argued, the requirement that the initial wave function changes 
slowly in the semiclassical limit is more difficult to justify. The main problem is that the 
Landau-type wave function (||D has the term B cos8y 2 /2h in the exponent and both terms, 
one coming from the initial wave function and one from the semiclassical Green function, 
have the same dependance on %. In this case one can proceed as follows. It is easy to check 
(see fln|)) that the requirement of the smoothness of the initial wave function is equivalent 
to the condition 

(3 < d 2 S(y, y')/dy 2 = m n /m 12} (13) 
where (3 = B cos 6 and the same condition for the derivative over y' . 
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But the value of the magnetic field in dimensionless units is very small. On the other 
hand for generic chaotic systems all elements of the monodromy matrix should be of the 
same order and there is no general reason why the above ratio should be small. Therefore 
we expect that for hyperbolic periodic orbits with large Lyapunov exponent the condition 
(|T3|) will be satisfied. We have checked numerically that the ratio mii/(mi 2 /3) grows quickly 
with increasing m\\ though for stable and almost stable orbits it can be of the order of 1. 

To treat the latter case a slightly different approximation can be used. Let us expand the 
action S(y, y') up to the second order and perform the integration over y and y' in (||) taking 
into account explicitly the dependance of the initial wave function on these variables. 
After simple algebra we get the same expression as in (|12|) but with the pre-factor D(r) and 
the action S(r) = rS p + AS(r) substituted by 

D(r) = m 2 i(r) + i(3(m n (r) + m 22 {r)) - f3 2 m 12 (r), 
A *S(r) = -^—r{yl(m 22 {r) + i^m 12 {r)) 

+ ^ 2 (m n (r) + i(3m 12 {r)) + 2y y' ). (14) 

These formulae are valid provided the third and higher terms inn the expansion of the action 
are small in comparison with quadratic terms. 

When j3 — > this result leads to Eq. fll^) . But for stable and almost stable regions where 
condition ( Dp is not well satisfied Eqs. (^) with ([14]) continue to be a good approximation 



for T in Eq. (|9]). In such cases one can (exactly or approximately) split the Hamiltonian 
into the sum of two Hamiltonians for longitudinal and transverse motion and the energy 
of the former is much bigger than the latter. Under these conditions (i) the period of 
motion is defined mainly by the longitudinal motion and (ii) the transverse motion is close 
to the motion in a quadratic potential similar to (|I]) with 9 = for a fixed time. But 
it is known |16j that the semiclassical time-dependent Green function for the motion in 
a quadratic potential coincides with the exact Green function and both are proportional 
to the exponential of a function quadratic in its coordinates. This means that for stable 
and almost stable regions there exist good reasons why higher than quadratic terms in the 
expansion of the action on transverse coordinates should be small and, therefore, Eqs. ( |T2"D 
and flT4]) are good approximations to the triple integral in (H). Of course, near bifurcations 
when by, definition, the quadratic form in the exponent is degenerate these equations require 
modification. 

There is an interesting limit of the above formulae. It corresponds to very clean devices 
where To is small and the probability of the tunneling through the second wall is much bigger 
than through the fist one. In this case R equals the reflection coefficient form the second 
wall and using the relation |r| 2 = 1 — \t\ 2 one concludes that the factor |t 2 *°^| 2 /(l — |-R| 2 ) 
in Eq. ([12]) tends to 1. One can show that all formulae in such a limit coincide with the 



semiclassical limit of a simplified model of resonant tunneling discussed in |TTJ and [[12] which 
is very convenient for numerical simulation. This model is based on the Bardeen transfer 
matrix |13[ according to which the probability of tunneling (or the imaginary part of the 
energy level) is given by 

r = 2vr]T|wg 2 <K£-£n), (15) 
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where E n denotes the exact energy levels in the quantum well and the coefficients W n are 
the matrix elements of the current between the wave function in the first well (ipi) an d the 
exact wave function (i/; n (y)) in the QW 

with the z-component of point q = (y, z) somewhere inside the barrier. 



Rewriting Eqs. (p~5f ) and (pT6[ ) through the Green function and its z-derivatives and using 
the same arguments as above it is possible to demonstrate that the semiclassical approxi- 
mation to the resulting expression are given by Eqs. (0) with \t { * ot) \ 2 /{l - \R\ 2 ) = 1. 

The same formulae can be used to derive an approximation similar to the Miller modifi- 



cation of the Gutzwiller trace formula |17] for a stable orbit. In such a case the monodromy 
matrix elements for r repetitions around a primitive periodic orbit are functions of exp(iur) 
where exp(±iu;) are eigenvalues of monodromy matrix. By representing Eqs. ( ]T2| ) and ([14]) 
as a power series J2 m exp(— iurm)C m and performing the summation over r one gets an 
expression of T as a sum over m: 

_ 7r(/?p)V* ri p_ , 2 2 

i_ l + (3p T 6XPl 2(l + (3p) [yo + y ° ))lVpl 

?24(|ri) m ™^^ 

J2S(E-E N>m ), (17) 

JV 



X 



X 



where p = m 12 / sincu, A = ft\J~pj (f3 2 p 2 — 1) and £7jv, m are approximate eigenvalues defined 



as m 



T7J: S p (E N)m - e±(m)) = 2nh(N + p p /A) and e ± = ^(m + 1/2). 
The index m has the meaning of the number of states in an effective transversal Hamil- 
tonian 

H(p,y) =uj(pp 2 + y 2 /p)/2. 

and the resulting expression corresponds to Eqs. (fl5|) and (^6[) in which one approximates 
ip n (y,z) by a product ip(z)ip m (y) where ip m (y) are eigenfunctions of this Hamiltonian. By 
truncating this sum at a finite value of m one obtains the Miller-type approximation which 
generalizes the tori quantization discussed in 



In conclusion we have developed a simple semiclassical theory of resonant tunneling. In 
the strict semiclassical limit the tunneling current is due to only special periodic orbits in 
the QW which are perpendicular to the LH wall and hit the RH wall of the QW, and their 
contributions are proportional to (m 2 i)~ 1 ^ 2 . For almost stable orbits one has to use Eq. ([14]) 
but the period of oscillations will deviate slightly from that of periodic orbits. For big stable 
regions the Miller-type formula ( |i~7D is appropriate. 
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